Remarks on astrophysical origin of the knee in cosmic ray spectrum 
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All astrophysical explanations of the knee in cosmic ray spectrum accept the hypothesis of its existence a priori 
without any doubt. But, there exist experimental evidences against this hypothesis. Experimental data on the 
knee in various secondary cosmic ray components obtained by KASCADE array as well as by many other arrays 
reveal the existence of the knee (A7 ~ 0.4 ■£■ 0.5) in electromagnetic component (e and 7 including Cherenkov 
light) and the absence of it or a very small knee (A7 m 0.1 -f- 0.2) in muonic and hadronic components. But, one 
could expect the inverse relationship in case of the astrophysical knee origin. A brief review of experimental data 
concerning this problem is given. 



1. Introduction 

The astrophysical origin of the so-called " knee" 
in primary cosmic ray spectrum in PeV region is 
very popular and is widely exploited in practice 
for a long time of about 50 years. This hypothe- 
sis seems to have begun a new life in the recent 
years after discovering some new Supernova rem- 
nants. But, experimental data accumulated up 
to now in cosmic ray physics seed some doubt on 
the "knee" existence. One could argue that there 
exists no theory able to explain pure power law 
spectrum in a very wide range from 10 10 to 10 20 
eV. But, such a theory does exist pp. It was pub- 
lished many years ago and proposed a universal 
mechanism of cosmic ray acceleration in the cos- 
mic plasma pinches. This theory can explain not 
only power law spectrum but even the value of 
its power law index 7. This mechanism can gen- 
erate cosmic rays up to the highest energy with 
the integral exponent 7 = \/3 « 1.73. Another 
outgoing feature of this mechanism is as follows: 
the index does not depend on a particle mass or 
charge. It is really a universal value. 

2. Brief review of the experimental data 

Let us look at the problem from the experi- 
mental point of view [2]. Direct measurements 
do not show any significant change in the spec- 



trum slope up to energy 1 PeV. The "knee" in 
the size spectrum is only seen in secondary cos- 
mic ray components measured by the indirect 
method of Extensive Air Shower (EAS). If pri- 
mary spectrum index is 7 and a secondary x- 
component depends on primary energy Eq as 
N x {E a ) ~ Eq then the EAS-size distribution on 
N x is P(N X ) ~ N~P, where (3 = 7/a. If the 
"knee" in primary spectrum does exist (let it be 
A7 = 0.5) one can predict with rather high accu- 
racy a relationship between the "knees" in all de- 
tectable secondary components: electromagnetic, 
muonic and hadronic. Index a commonly used 
for electrons is equal to a e — 1.15 — 1.25 (so 
A/? e = Aj/a e w 0.42) and that for muons and 
hadrons is equal to a^.h — 0.8 — 0.95. 

Therefore, the visible "knee" w 
0.5/0.87 « 0.57) in muonic and hadronic compo- 
nents is expected to be bigger than in electrons (in 
EAS size spectrum) and even bigger than in pri- 
mary spectrum, because a Mi h < 1 while a e > 1. 
It is a very strict limitation because a^^, and a e 
were measured in many experiments and were cal- 
culated by many authors using different Monte- 
Carlo programs. In this paper I would like to 
emphasize that experimental data on a visible 
"knee" in different secondary components con- 
tradict this limitation. In other words, the val- 
ues of the "knee" in primary cosmic ray spec- 
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Table 1 

Measurements, predictions or expectations for the knee value 
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trum reconstructed from different secondary com- 
ponents are not equal and contradict each other 
as is seen in Tabic 1. In all experiments shown 
above A/3 e > ps A@h- The most detailed 

and careful measurements in the knee region have 
been performed by KASCADE experiment spe- 
cially designed to solve this problem. The "knees" 
claimed by the authors in muonic and hadronic 
components are much less than expected ones 
(even equal to for 'e-poor' EAS's) and are very 
smooth and thus could be regarded from my point 
of view as a smooth change of the power index 
in a very wide range due to methodical or other 
reasons to be discussed below. Therefore, their 
values could be regarded as the upper limits for 
the "knee" . If so, the KASCADE data are consis- 
tent with the "zero knee" hypothesis made on the 
basis of all available data analysis [5] of the deep 
underground muon experiments (Baksan, Frejus, 
NUSEX, MACRO, SOUDAN2). In fact, only 2 
experimental groups have confirmed the existence 
of the knee in the EAS muonic component compa- 
rable with that for the electron component: MSU 
EAS and EAS-TOP. But their data need special 



comments. MSU EAS muon data obtained with 
muon detectors do not reveal a knee. What is 
usually shown by the authors is a muon number 
spectrum recalculated from the electron number 
spectrum. Therefore, it must have the knee close 
to that in EAS size spectrum. As for EAS-TOP 
data, their muon number spectra shown in 8 
(and in other publications as well) could be suc- 
cessfully fitted by a pure power law by taking into 
account statistical errors and the points spread. 
Moreover, the results of these two experiments 
contradict each other because EAS-TOP claimed 
the knee in muonic component (E^ > 1 GeV) at 
7V M = 10 4 - 9 while MSU EAS group (E^ > lOGeV) 
do not see any knee above = 10 4 (the differ- 
ence in muon threshold energy gives a factor of 
5-6 while the difference in figures is higher and 
the difference in altitude of ~ 200g/cm 2 is neg- 
ligible for muons). The result of AKENO com- 
bined in [Hj with SUGAR data JI2] gave a rather 
strict limitation for muon size spectrum as fol- 
lows: " ...there is no change of slope in muon 
size spectrum between 10 5 and 10 9 ". Another 
strict limitation starting from a very low muon 
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size has been given by GAMMA data 13 for 
mucins with energy above 5 GeV: "The spectrum 
above Nn > 7000 is well described by the power 
law dependence F(> N^) ~ N~ 2 - 18 ". 

An absolute value of the power law index (fi^) 
has to be discussed separately As it was demon- 
strated in |17j the indices of primary spectrum 
derived from measurements by different experi- 
mental groups are in disagreement. In our ap- 
proach 7 = const and everything is explained by 
the index a behaviour, which is not constant in 
different experiments. This can be explained at 
least for the measurements of primary spectrum 
through the muonic component. As it has been 
shown |18j muon number spectrum is connected 
with the primary cosmic ray spectrum and fol- 
lows power law if primary spectrum does so. In 
case of an infinite area detector, the simplest one 
corresponding to the total number of muons in 
EAS when a M w 0.7 0.8 (let it be equal to 
0.7 for estimations), the index of muon number 
distribution in asymptotics should be equal to 
(3fj, = j/ctfj, sw 1.7/0.7 w 2.4. Low multiplicities 
(m < 5) are produced by low energy primaries 
and can be described by Poissonian or binomial 
distributions characteristic for the discrete values. 
Such a spectrum was observed in many under- 
ground experiments (see for example [H]) using 
large area detectors for recording of muons with 
a high energy threshold without any EAS trigger 
(case (a)). A similar result is expected for small 
underground detectors but only the ratio of sin- 
gle muon flux (m=l) to that for muons in groups 
(m > 1) is expected to be much higher. Another 
interesting feature of such experiments (without 
EAS core location) is the so-called core attraction 
effect [20]: the higher multiplicity in the detector 
is - the sharper recorded EAS core lateral dis- 
tribution centered in the detector is. This effect 
results in a measurement of not the total number 
of muons integrated over all distances up to infin- 
ity but of only the central muon density having 
a pole in the axis where the exponent is equal to 
otfx i=i 1. Therefore, the measured muon num- 
ber spectrum at very high multiplicity (where 
the pole does exist) should be flatter [2] hav- 
ing the index ^ = 7/0^ w 1.7/1.0 = 1.7. The 
latter could probably explain the problem with 



the highest multiplicity events observed in the 
CERN detectors ALEPH and DELPHI. In real 
underground experiments due to different experi- 
mental and trigger conditions and due to different 
statistics, the measured integral index should lie 
between these two extreme values {j3 ~ 2.4-^ 1.7). 

The opposite case (b) occurs when the muon 
detector is working under the EAS array trigger. 
If the data are corrected for the trigger bias then 
the EAS muon size spectrum will follow the power 
law with j3 « 2.4. If no corrections are applied 
then the measured muon multiplicity spectrum is 
very flat at low multiplicity w 0.7) |18ll9l20j 
according to the Poissonian distribution with the 
high mathematical expectation (mo > 1) which 
depends on the EAS size threshold and on the 
muon detector total area. Only at a high enough 
number of muons, the spectrum becomes as steep 
as in case (a). This effect could imitate the ex- 
pected knee in primary spectrum. If muon den- 
sity is measured at fixed core distances (as in 
KASCADE and many other EAS experiments) 
then the core attraction effect does not work and 
no flattening occurs at high multiplicity in such 
experiments. Here could be probably found the 
difference in experimental data at the highest 
muon density measurements between EAS arrays 
and underground detectors. 

Unfortunately there exist very few data on the 
hadronic EAS size spectrum measurements free 
of trigger conditions bias. I found only one ex- 
periment (0]) providing such data. 

We can conclude that none of the experiments 
has shown undoubtedly the knee existence in 
muonic or hadronic component in the needed 
place and of the needed size as could be expected 
in case of real knee existence in primary spec- 
trum. 1 

3. Summary 

Summarizing the experimental data shown 
above we can conclude that the question "Does 
the knee exist or not?" is still open. Therefore, 
the main efforts should be directed now not to 
explanations of the "knee" but to the question 

1 concerning the KASCADE muonic data this has been 
first mentioned in 1221 . 
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above. Only when this question is answered posi- 
tively an astrophysical explanation of the "knee" 
existence in primary cosmic ray spectrum will be 
needed. Up to date I incline to a negative an- 
swer to this question, because only in this case 
we could understand the experimental data. In 
the frame of the approach proposed in [2] there 
exists a natural explanation of the visible "dis- 
agreement" of the experimental data shown in 
Table 1: the "knee" does not exist in primary 
cosmic ray spectrum and should be observed only 
in electromagnetic component due to non-power 
law behaviour of the N e (Eo) in a region of ~ 100 
TeV/nucleon when hadron core reaches the sea 
level and coreless EAS's become coreful or normal 
ones. Muons and hadrons are almost insensitive 
to this critical point of the EAS development and 
follow a power law much better. 

Looking at Table 1 we can conclude that the 
measured data are more consistent with the pre- 
diction [2] shown in the last line than with the 
expectations made above under a supposition of 
the "knee" existence (first line). And finally: if 
the answer is negative then we probably have a 
natural universal theory of cosmic ray accelera- 
tion [1]. To confirm or to reject this theory we 
have to measure the spectrum index 7 with high 
accuracy and answer the question: " Does it equal 
to V3 or not?" 
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